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ABSTRACT 

A 8«t of rate equations including strong turbulenoe effects and anomalous 
resistivity are solved using parameters which model several solar type III 
bursts. The electron exciters of these bursts have been detected at earth 
orbit. The analysis has enabled us to provide quantitative comparisons 
between several of the observed phenomena and the theory. Using an analytic 
model for the time evolution of the energetic electron exciter, we find that 
the exciter distributions observed at 1 AU are unstable to the excitation of 
the linear bump> in-tail instability, amplifying Langmuir waves above the 
threshold for the oscillating two stream instability (OTSI). The OTSl, and 
the attendant anomalous resistivity produce a rapid spectral transfer of 
Langmuir waves to short wavelengths, out of resonance with the electron 
exciter. Further energy loss of the beam is thus precluded. The various 

parameters needed to model the bursts are extrapolated inside 1 AU with 
similar results. Again, the OTSI is excited and decouples the electron beam 
from the Langmuir radiation. Reabsorption of the Langmuir waves by the beam 
is shown to be unimportant in all cases, even at 0.1 AU. 

The theory provides a natural explanation for the observed relationship 
between radio flux, 1, and the electron flux, Jg. When the OTSI is weakly 
excited the theory predicts that I«Jg» as observed; while for stronger bursts, 
the theoretical result is with o»2.i1. We find that the value of a is 

only weakly dependent on variations in beam density and the spectral index of 
the exciter distribution function. In all cases the theoretical results agree 
closely with the observed values. 


I. MIRODUCTIQM 


Observations of solar type III bursts have reached a level of sophlstl- 
oation which permits both qualitative and quantitative comparison between 
theory and experiment. In the accompanying paper (Smith, Goldstein and 
Papadopoulos 1979, henceforth referred to as Paper III— Papers I and II being, 
respectively, Papadopoulos, Goldstein and Smith 1974 and Smith, Goldstein and 
Papadopoulos 1976), we presented a theory of strong plamna turbulence as it 
applies to type III burst phenomena. In the present paper we extend the 
analysis of Paper III in am effort to provide quantitative comparisons between 
the strong turbulence theory and many aspects of type III phenomenology. 

In §11, the basic equations of the theory are solved numerically for 
parameters appropriate to several specific bursts. The solutions demonstrate 
the importance of the role played by the oscillating two stream instability 
(OTSI) and anomalous resistivity in decoupling the electron exciter from the 
plasma turbulence. Reabsorption, included in the analysis, is shown to be 
unimportant in all calculations at all heliocentric distances. 

The computed levels of Langmuir turbulence agree well with the recent 
reports by Gurnett and Andersen (1976 and 1977) of intense sporadic bursts of 
electrostatic noise observed in association with type III bursts. The 
predicted intensity of electromagnetic noise resulting from the Langmuir 
turbulence is found to be comparable to the observed levels. In addition, we 
show that the puzzling correlation between electron flux and radio intensity 
(Fitzenreiter, Evans and Lin 1976) has a natural explanation in the functional 
dependence of the OTSI growth rate on levels of Langmuir energy density. 
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II. miMERICAL CALCaLATlOMS 


The model rate equations describing the spectral evolution of electron 
plasma oscillations, including the effects of the linear beam-plamaa instab- 
ility, the oscillating two stream parametric instability and anomalous 
resistivity were first derived in Paper II and were more fully described in 
Paper III. The integration of the system of equations (III-3«44)-(IIl-3>47) 
is the subject of this section. [Roman numerals refer to the paper in which 
the equations are found.] Definitions of the various growth and damping rates 
(Yl» Yqjs* "^hL ^ found in Paper III. The integration 
required discretizing wavenumber space, or equivalently - phase velocity. 
Thus, there were of the order of 200 coupled non-linear differential equations 
integrated in time over the interval -0.7 ^ w/kc ^ 0 . 1 . 

The starting condition (tsO) was the arrival at the location of Interest 
(a point between 0.1 and 1 AU) of energetic electrons with streaming veloci- 
ties centered near 0.7c. The exact velocity distribution of the electron beam 
at that location was given by the beam evolution model described in SIV of 
Paper III. That beam evolution model was first constructed for a single 
event; it has but one free parameter, L — the path length traveled by the 
exciter. By varying L, the time behavior of the model at 1 AU could be 
adjusted to fit observations by Lin and coworkers of several different bursts 
(Lin 1974, and Fitzenreiter , Evans and Lin 1976). In all cases, the observers 
found that at first only the fastest particles were detected (with energies 
</*100 keV). This was followed over an extended period of time (>^20-30 minutes) 
by the arrival of an increasing flux of slower particles (cf. Figure 1). The 
long time over which the electron flux was observed to increase indicates that 
the particles experienced significant amounts of pitch-angle scattering; more 
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than can b« aooountad for by interactlona with Langmuir turbulanoe axoitad by 
tha bump» install Instability. In oonstruoting our modal of tha baam avolution 
(of. Papar III), wa took this into aooount by raquirlng that our phanomano- 
logieal modal raproduoa tha tima bahavlor of tha alaotron flux that Lin and 
ooHorkars obsarvad. ^a physical origin of this enhancad pitch-angle scatter- 
ing is probably due to the interaction of the electrons with the magnetic 
fluctuations known to permeate the Interplanetary medium. In this context, it 
is well to point out that calculations involving solutions to the time 
dependent quasilinear (i.e., weak turbulence) equations have been published by 
Takakura and Shibahashi (1976), Takakura (1977), and Magelssen and Snith 
(1977) in which the effects of magnetic pitch-angle scattering were ignored. 

In Magelssen and Smith's calculation, which attempted to model the type 
III burst of 16 May 1971, the calculated electron distribution evolved from 
100 keV to below 20 keV in less than 5 minutes (v. Figure 13 of their paper). 
As a consequence of this rapid evolution of the exciter distribution, reabsorp- 
tion of the linerly excited Langmuir waves was found to be very important. In 
a situation in which the flux of low energy electrons is rapidly increasing 
with time, reabsorption occurs as the peak of the electron distribution 
function, initially located at rather large velocities, shifts to lower 
velocities. Langmuir waves that had been initially excited when the slope of 
the electron distribution function was positive soon find themselves inter- 
acting with the portion of the distribution with negative slope. The waves 
are then quickly damped or reabsorbed. 

In our analysis reabsorption will be unimportant as long as the beam 
evolves slowly compared to the time necessary for the OTSI to effect a 
spectral transfer of the Langmuir radiation to large wavenumbers, out of 
resonance with the exciter distribution. In situations where the beam 
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evolution is as rapid as that considered by Magelssen and Smith (1377), 
reabsorption can be an important eflcct, however, as shown in Figure 1, Lin's 
observations indicate that more than 30 minutes actually elapsed before the 
peak of the distribution evolved to 30 keV. Consequently, Magelssen and 
Smith's calculation greatly overemphasizes the effects of reabsorption. As we 
show below, the OTSI produces a spectral transfer of the Langmuir radiation on 
time scales of less than a second, and thus reabsorption is unimportant in our 
analysis . 

As an initial condition in our calculations, the Langmuir and ion waves 
are taken to be in thermal equilibrium with the solar wind plasma, the 
equilibrium energy density of Langmuir waves being the enhanced noise levels 
produced by the suprathermal tails of the electron distribution function (v. 
eq. III-3.41). 

The rate equations were integrated using a predictor-corrector algorithm 
with variable step-size (Eiserike and Silver 1971). For each of several 
bursts the code was run at heliocentric distances ranging from 0.1 to 1.0 AU. 
In the following discussion, the results of the calculations are described for 
each burst in turn; first for parameters vdiich model observations at 1 AU and 
then for parameters which model observations between 0.1 and 1 AU. 

(a) The burst of 16 Mav 1971 at 1 AU. 

-17 -2 

This was an isolated and fairly intense burst (reaching nearly 10 W m 
Hz ^ — Fitzenreiter , Evans and Lin 1976) observed by instruments on both the 
IMP-6 spacecraft and the Apollo 16 subsatellite. These two spacecraft were 
able to simultaneously measure both the radio and particle fluxes. [Electron 
spectra, estimates of the beam densities, and estimates of the path lengths 
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tr«v«rs«d by the beta before reeohing 1 AU for this and the other bursts were 
kindly provided to us by R. P. Lin snd R. J. Fitzenreiter (slso see Lin 197A).3 

The local plaass frequency, f , at 1 AU on May 16, 1971 was about 30 kHz. 

pe 

Electrons with energies in excess of 100 keV and radio noise at 55 kHz (*2f_.) 
were first detected by instruments on lMP-6 at 1305 UT (see Figure 1, vdiich is 
taken fr«n Lin et al. 1973) • By 1335 UT the radio noise had reached its 
maxiauo and little further evolution of the electron spectrum was observed. 
Between 1305 and 1312 UT the peak of the electron spectrum had shifted from 
several hundred keV to about 80 keV. Subsequent evolution was slow, and it 
was not until 1324 UT that the peak had moved to about 50 keV. By 1332 UT the 
peak was at 40 keV, and had evolved to 30 keV by 1335 UT. At later times no 
clearly defined peak was present. Thus, thirty minutes elapsed between the 
time that 100 keV electrons were first observed and the time that the peak of 
the spectrum had moved down to 30 keV. As discussed above, our model of the 
beam is constructed to have this slow evolution. 

For the burst on 16 May 1971 (at 1 AU), the path length L was estimated to 
be 1.5 AU. The ratio of the density of the beam to that of the solar wind was 
n“5x10”^, and the spectral index of the electron distribution was c»4.6. 
These parameters are summarized in Table 1 for this and all other calcula> 
tlons. Below the peak velocity, Bp(t), the distribution function, f(e) had a 
positive slope throughout the 30 minutes between 1305-1335 UT. With these 
parameters, the rate equations were integrated and the evolution of the plasma 
Instabilities followed in some detail. 

The solution is presented graphically in Figure 2, where f,j.(6), W(+k) , 
W(-k), and (6n/n)*=S(k)(kX )* are plotted on a logarithmic scale against 
^ph*“k^*^® at various times (Fig. 2a-2f). 
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Initially the linearly unstable beam produces resonant plasma waves 

(Indicated by cross hatohlng In Pig. 2a) that grow until the turbulence level 

exceeds the threshold of the OTSI (eq. III-3.<:9) At this time W*>1x10*^ 

(corresponding to an electric field of 2 mV/m). With the onset of the OTSI, 

aperiodic Ion waves are excited (the gray shading In Pig. 2c>2f) at the growth 

rate given by equation (III-3.28) This Is also the rate at which the "pump" 

waves earlier excited by the linear Instability are scattered to shorter 

wavelength "daughter" waves propagating parallel and antiparallel to the beam 

(Fig. 2b-2d). If there were no other physical processes operating, the pump 

-5 

waves would saturate near the thresholu value of W«10 . Although this 
saturation level is well below that expected on the basis of weak turbulence 
theory (W^/^l), resonant energy exchange between the beam and the pump waves 
would continue, eventually leading to substantial deceleration of the beam. 

However, the OTSI does not stabilize at this enhanced level; for as the 
ion turbulence grows, the ion density fluctuations modify the Bohm-Gross 
dispersion relation so that the threshold of the OTSI is lowered (cf. eq. 
111-3.3^). This further reduces the energy in pump waves (Fig. 2e) . Finally, 
the a-c anomalous resistivity associated with the correlated ion fluctuations 
scatters all long wavelength Langmuir waves to short wavelengths (Fig. 2f ) . 
In our calculations the collapse to short wavelengths ceases when Landau 
damping by the thermal electrons balances the scattering. At this point the 
process has stabilized. The electron beam is now .completely decoupled from 
the plasma waves, and no further energy exchange will take place between the 
beam and the radiation field. If we were to continue the calculation, which 
we have not done for reasons of economy, we would expect that gradually the 
ion fluctuations and Langmuir waves would simultaneously decay. This is a 
process that has not been fully investigated theoretically, but should take 


8 



pUoe on th« time scale of Ion Landau damping (eq. ZIl-3.39). Eventually, 
following tbe deoay of the ion fluotuations , the linear instability oould 
begin again. Suoh a oyolio process, similar to relaxation oscillations, would 
continue until the electron beam merges wii?i the ambient sol'ar wind electron 
distribution, at tdiich time the total distribution function f^(0) would not 
have a positive slope. [He reported similar behavior in an earlier version of 
these calculations (Paper II). Recently, H. Rowland has observed the same 
cyclic phenomenon using Vlasov type simulations for both electrons and ions 
(Rowland, private communic. ^.ion) . ] Merging of the exciter distribution 
function with the solar wind was observed by Gurnett and Frank (1975) during 
the type III burst on 5 November 1974. Another mechanism for ending the 
instability would be the disappearance of the velocity space anisotropy. In 
fact, this latter possibility seems more consistent with observations of the 
anisotropy made by Lin during the 16 May event (Lin, private communication). 

Little more than 0.1s elapsed between the time of OTSI threshold was 
reached and the final stabilization of the various instabilities. In 0.1s the 
electron distribution is essentially constant, vindicating our claim that 
neither reabsorption nor quasillnear relaxation is Important. Even the 2.8 
minutes during which the linear instability operates at 1 AU is insufficient 
time for the beam to reabsorb an appreciable amount of wave energy. We 
reiterate that our calculation Includes the damping of the waves (or reab- 
sorption) at all phase velocities for which 3f/3B<0 (cf. eq. III-3.43). 

(b) The burst of 28 February 1972 at 1 AU. 

This burst was less intense than the one on 16 May 1971, the electron flux 
was lower, the distribution function had a softer spectum (;=10.4), but. 
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nonetbelessi the beam waa relatively dense (riBl.4x10 ). The plasma frequency 
at 1 AU was .^30 KHz, and the ambient electron temperature was >1.2x10^*K. The 
beam is estimated to have traversed a path length L«1.9 AU. At 1 AU stabili- 
zation proceeds in much the same way it did for the 16 May 1971 event (Fig. 
3a-3o). However, because c is larger, threshold for the OTSI occurs at 

smaller 6-0.31 with W-5x10“®. 

P 

As before, stabilization proceeds rapidly once OTSI threshold is attained, 
and neither reabsorption nor quasilinear relaxation are important. Figure 3a 
shows the pump waves near threshold; in F. ^’ure 3b the OTSI has scattered most 
of those waves to lower phase velocities .both parallel and anti-parallel to 
the electron beam); and in Figure 3c stabilization is complete. Note that the 
electron flux at stabilization is significantly smaller for this and the 25 
May 1972 burst, discussed below, than it was for the burst on 16 May 1971 
(Fig. 2). 

(c) The burst of 25 Mav 1Q72. 

This was similar In some respects to the one on 28 Fe: -uary 1972 Just 
discussed. The spectrum was again soft (;s6.4), the electron and radio fluxes 
were rather modest compared to the 16 May 1971 event, and the electron beam 
was even less dense (n-5.6x10”^) than on 28 February 1972. The ambient 

electron temperature was 1.2x10**K, Lsl.5 AU, and f s22 KHz. At 1 AU the 

p6 

-5 

OTSI threshold was reached at Bp-0.33 and W-4x10 . Stabilization again 

proceeded quickly (Fig. 3d-3f) . 

(d) Variation with heliocentric distance. 

Our discussion has thus far been confined to calculations using parameters 
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ohariottrlstio of oondltions ot 1 AU, booaus* it it only tt 1 AO that siaul- 
tanaoua radio and partiole data hava thus far baan availabla. Obsarvationa of 
plaaaa and radio wavaa, but not sinultanaoua alaotron BaaBuramanta, hava baan 
publisbad following tha HBLIOS 1 and 2 Bissions to 0.A5 AO. This has proBptad 
us to invastigata solutions to tha rata aquations at various haliooantrio 
distanoas batwaan 0.1 and 1 AO. 

At 0.5 AO, for the 16 May 1971 event, we took the plasoa frequency to be 

60 KHz, U0.7 AO, and assisned that c and c were constant with heliooentrio 

distance— (Table 1). The results are qualitatively siailar to those obtained 

-5 

at 1 AO, although the threshold energy density is now H>3x10 (corresponding 
to an electric field of bV/b). Because one is closer to the corona, the 
electron distribution evolves more rapidly so that the peak of the spectrua is 
near 0.M5c tdien the OTSI threuhold is reached. The more rapid beam evolution 
at smaller heliocentric distances implies that reabsorptlon of the electron 
oscillations by the beam probably plays a somewhat more important role in the 
inner solar system than at 1 AU. Monethless, as shotm in Figure Ma-4o, 
stabilization proceeds via the OTSI and anomalous resistivity, and not by 
reabsorption. [Note that in Figures 4 and 5, for reason of economy, the 
calculation have not been carried through to the point where the OTSI is 
completely stabilizaed by anomalous resistivity (cf. Figs. 2f, 3c and 3f).l 
Figures 4d»4e show the results of integrating the rate equations for the 
16 May 1971 burst at 0.1 AU. For this calculation f_.>!300 KHz and LsO.1 AU. 

pe 

stabilization proceeds as before via the OTSI, and again reabsorption is 

.q 

unimportant. At threshold, H*1»10 and S^*0.3c. 

Similar calculations were performed for the 25 May 1972 and 28 February 
1972 events at both 0.5 and 0.1 AU. These results are shown in Figure 5 (0.5 
AU) and Figure 6 (0.1 AU). Again, these bursts reached the OTSI threshold at 
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■oa«Mh«t lower olootron volooltios than ms the o«s« with the sore intense 
burst of 16 Itaf 1971 (of. Fig. 4). Nonetheless, even st 0.1 AU reabsorption 
is not signifiosnt, and the bees beooses deooupled froa the Lsngauir turbu* 
lenoe after snotislous resistivity stabilizes the OTSI. 

These three bursts (16 May 1971, 25 May 1972 and 28 February 1972) are the 
only ones for tdiioh we had eleotron data sufficiently ooaplete to provide the 
paraaeters needed in the beaa evolution aodel, viz. L, n, and c. They are 
included in the events analyzed by Fitzenreiter, Evans and Lin (1976). In 
that analysis, Fitzenreiter, et al. discovered a relationship between the 
radio flux, I, and eleotron flux, Vhen Jg was less than 100 

(oa**8»ster)“\ I and Jg were related by I(2ui^)«Jg® with o*1. For Jg>100 
(cB*»s*ster)“\ I(2w^)«Jg® with o«2.t. 

The three bursts we were able to aodel represent examples in which was 
either of order 2.4 or changed from 1 to 2.4 during the event. The 16 May 
1971 burst was an exaaple for which o«2.6 frca the onset of the burst at 1 AU 
(1305 UT). For both the 28 February and 24 May 1972 events, Jg was less than 
100 (ca**s*ster)~^ for soae ten minutes after (wiset, and a*1. Subsequently, 
Jg exceeded 100 (cm**s*ster)~^ and a increneed to about 2.4. He will return 
to discuss these two classes of events in more detail following the discussion 
of solitons, which follows. 


111. SOLITOMS 

In ordtr to giro • tbtorotioal dosoription of th« •Itotrooagnetio rodia* 
tion propoFtloa of thaaa bursts and tbs oorralatlon batwaan and 1, wa oust 
oonsldar tha ^atial struotura of tha oollapaad short wvalangth plasma 
osoillationa raprasantad in Figura 2f. 

A comparison of Figures 2a and 2f illustrates tha collapse of the long 
wavelength pump wave into short wavelength daughter wavM. As pointed out by 
Manhaimar and Papadopoulos (1975)t the aquations that describe the OTSI are 
tha Fourier representations of those that describe soliton formation. In 
configuration space tha pump waves comprise a wavepaeket with dimensions 
AZBl/Ak^«1200x^a8 km with x^«700 on at 1 AU. Following collapse, the scale 
length of the solitons, Ax , is Ax >50x km and W has reached nearly 

9 8 9 

10 . This is something of an overestimate oecause our numerical calculation 

does not take into account the effect of the magnetic field in inhibiting 
collapse in the transverse directions. When U<<1 each wavepaeket collapses 
into one soliton, so the spacing between solitons is approximately az. 

Soliton formation oan be investigated directly from equations (In-3.11)- 
(III-3.12) (Zakharov 1972, Budakov 1973) • Papadopoulos and Frejnd (1978) have 
argued that the solitons collapse only in the direction parallel to fi, in 
which case the energy density in the soliton can be estimated to be (Kingsep, 
Budakov, and Sudan 1973* Degtyarov, Zakharov, and Budakov 1976). 

V • 12(kX^)* (1 T,/Tj) - 24(kx^)* 

With kX^-SwIO”^ (cf. Fig. 2f), one finds W>10"^, which corresponds to an 
electric field on 60 nV/rn. 
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Until recently, there has been little rxperiaental evidence that such 
collapsed, intense electric fields aoooapany type III bursts. The reason is 
clear. In a 400 ka/s solar wind, a 350m Mliton is convected past a 30a 

dipole antenna in little acre than a Billisecond. Plasma wave experiments on 
spacecraft have been designed with electronic response times no faster than 
20ms (D. A. Gurnett, private caoBunicatlon) so that the likelihood of detect- 
ing a fully collapsed soliton along its parallel dimension is small. 

Papa(k>poulos and Freund (1978) noted that because there is little or no 
collapse sdong the transverse direction, solitons might be detected as the 
transverse dimension, Lj^, is swept past an auitenna. Depending on the ratio 
u /O , this size might be connected with either the beam electron Larmor 
radius, or the transverse correlation length of the instability, * 

Vg/Yj^. For example, at 0.45 AU, km, so that the 

transverse size of a soliton might take longer than 60ms to convect past a 
spacecraft, and would then be detectable. While the two HELIOS spacecraft 
were near 0.45 AU, Gurnett emd Anderson (1976, 1977) did report observations 
of intense, sporadic bursts of plasma oscillations in association with four 
type III bursts. The intense plasma oscillations at u> (£•'‘15 mV/m) were 
observed only axter the electromagnetic radiation at 2u^ had reached its 
maximum and started to decay. 

Because the strong bursts of plasma oscillations often lasted for 1/3s of 
longer, Gurnett and Anderson could not have observed a fully collapsed soliton 
convecting past the spacecraft as long as the magnetic field was approximately 
parallel to the flow direction of the solar wind. However, this observation 
is consistent with the transverse dimension, X, of a soliton being convected 
past the spacecraft if 25 km < X < 90 km, which satisfies the requirement that 
X be greater than (Papadopoulos and Freund 1978). When detailed magnetic 
field data from the HELIOS magnetometer experiments beccxnes available it 
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should be possible to determine the orientation of solitons when the spiky 
turbulence was observed. 

An alternative interpretation is possible, though without additional 

analysis it is at present only a suggestion, viz., that after the pump wave 

has been depleted, the collapse wave packet begins to spread diffusively at 

its group velocity. At 0.45 AU, X »400 cm and v =>10’ cm/s. In Is the packet 

e g 

would diffuse to a size of 7x10’ cm, which would take «r1/4s to convect past 
the antenna system. However, at these dimensions, expanding solitons would 
start interacting with each other. This situation has not been investigated 
and is beyond the scope of the present paper. 

(a) Radiation from Solitons 

A soliton is a spatially confined structure in which intense electrostatic 
fields oscillate at frequencies very close to As pointed out by Papado- 

poulos and Freund (1978), such fields will produce a second order current 
where n.| is the density perturbation consistent with Poisson's 
equation, and is the velocity perturbation produced by the soliton electric 
field. If that field has the form (in cylindrical coordinates) 

JE.(n»t) = e E(r,z) sin(o) t) 
z e 

then (Papadopoulos and Freund 1978) 

-e e /E 3£\ sin(2w t) 

si" \ / 

e 

Note that radiation will be produced directly at 2u)^. For the spatial 
structure of a soliton one can take (Degtyarov el al. 1976, Papadopoulos and 
Freund 1978) 
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B(r,i) • 8«oh (k^i) 

Pipadopoulos and Fraund (1978) then find that the total volugie ealaalvity, 
J(2«d^) la 


J(2d)g) 


3 /3yvj\7oE*y 1 y 
8 Ac / V8«A^^ ICqL/ 


(3.1) 


where, again, Az ia the parallel dimenaion of the linearly una table wave 
packet and is the wavelength of the elect r<»Dagnetic wave at 2u^. In 
writing equation (3.1) we have aaavaned that k^*L*»M, aa it la throughout the 
Interplanetary medium. The intensity of emlaaion Just outside a spherical 
shell of radius R and thickness AR centered on the sun is (Gurnett and Frank 
1975) 


I 


JR 

(2ui^/2it) 

6 


(3.2) 


Equations (3.1) and (3.2) can be evaluated as the rate equations are 
integrated. For the 16 May 1971 burst, at the time of soliton formation at 1 
AU (Fig. 2f), we find that Ij'IxIO”^'^ W m“^ s”^ with Az"3(v. /v )X «2 km and 

D G c 

Lb 24 km. This is close to the peak Intensity that was observed during that 
event at 55 KHz. As stated above, our calculation probably over estimates 
E*/8w in the soliton, and so, coupled with the uncertainties in determining 
the other parameters of the model, this seems to be quite satisfactory 
agreement . 


(b) Radio and Electron Flux Relationship. 

As we mentioned, Fitzenreiter, Evans and Lin (1976) discovered a strong 
correlation between simultaneously observed radio and electron fluxes. 
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Papadopoulos and Freund (1978) pointed out that this correlation results Vroa 
the scaling of W(k^) with Yqxs* parameter range of interest in the 

solar wind, this scaling law is (Papadopoulos 1975i Rowland 1977) 

^OTS ^ > “t 

W(k^) </> Ynxs '*^*'0^ ^ 

(3.3) 

The growth rate used in our computations (eq. 1II-.3.28) follows this scaling 

law very closely (Fig. 7). In fact, the asymptotic regime W(k^)>^YQ^g is 

reached for values of W(k^) somewhat less than m/M. 

o 

When the OTSI stabilizes the beam plasma instability, and \“”b' 

With I(2«^)«W(k^) (eq. 3.1 and 3.2), one has I(2u )«n. or n^. The electron 
flux, Jg, is in turn proportional to nj^<v>, vrtiere <v> is the velocity of the 
of the exciter distribution— it is that part of the electron beam that 
drives the linear instability. Thus, with f'/'v”^, one has 
so that 

I(2u ) « j (1-c)/(2-c) > W(k ) > W_ 

b Cl 0 1 

or 

I(2w ) « j 2(1-C)/(2 -c) j ^ 

e E 0 

(3.4) 

The 16 May 1971 event exhibited no transition and a was 2.63 during the 
entire event. We noted earlier that for this event c«4.6, so that, from 
equation (3.4) we calculate as2.77. 
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On 28 F«bruary 1972 the electron spectrum had a slope of c«10.4, which 
yields as2.24 late in the event. Fitzenreiter et al., observed acO.84 at 
Jg<100 (cm**a»ster)”^ anu as2.38 at larger Jg. 

Similarly, the 25 May 1972 burst had c<»6.4, implying theoretical values of 
as1.23 and 2.46, while asO.92 and 2.54 were observed. Not only are computed 
values of a remarkably close to those observed but we also find that the 
relative order is accounted for, i.e., a is largest on 16 May 1971 and 
asallest on 28 February 1972. These results are summarized in Figure 8. 

The scaling 4hat o should not be less than one. From 

Figure 7, however, it is clear that very close to threshold, W(k^)«rY^^g with 
v<1 . Thus, for bursts which initially only weakly excite the OTSI, one does 
expect values of o<1. In fact our numerical calculations confirm that for 
both of the 1972 bursts (28 February and 25 May), W(k^) does not exceed by 
nearly as much as it did during the 16 May 1971 burst. 
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IV. CQNCLUSIQtlS 


In this series of papers we have sought to demonstrate that the physics of 
solar type III radio bursts is a manifestation of strong pla^a turbulence 
processes. The oscillating two stream Instability together with anomalous 
resistivity provide the means by which electrons with energies from tens of 
keV up to the order of 100 keV are able to traverse distances greater that 
10' 'em with only modest energy losses. In modeling the nonlinear processes by 
means of the rate equations derived in Paper II we have been able to study 
type III bursts using parameters typical of those observed in the solar wind. 
This is often difficult to do in numerical simulations of parametric instabili- 
ties, primarily because of the very low beam densities found in the Interplane- 
-i| -6 

tary medium (n“10 -10 ). We have found that for the parameter range we were 

able to model, electron beams will produce levels of Langmuir turbulence which 
exceed the threshold for excitation of the OTSI. Furthermore, these levels 
are reached on time scales that are rapid compared to the time scale over 
which the peak of the electron distribution function of the stream evolves to 
lower velocities. Consequently, the spectral transfer of Langmuir radiation 
to large wavenumbers (small wavelengths), occurs rapidly. The phase velocity 
of the radiation then lies in a region of phase space in which there are few 
if any exciter electrons — the electron beam then propagates without further 
energy loses. Because of the relatively slow time evolution of the electron 
exciter, reabsorption of the Langmuir radiation was of minor importance. 

As pointed out by Papadopoulos and Freund (1978), theories based on weak 
turbulence predict far higher intensities of radio noise (by up to three 
orders of magnitude) than is observed. In contrast, because electromagnetic 
radiation from solitons is proportional to E*/8ir (rather than E'/(8'!r)* ) , the 
radiation levels computed are close to those observed. Furthermore, the 
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should ba possible to determine the orientation of solitons when the spiky 
turbulence was observed. 

An alternative interpretation is possible, though without additional 
analysis it is at present only a su^estlon, viz., that after the pump wave 
has been depleted, the collapse wave packet begins to spread diffusively at 
its group velocity. At 0.45 AU, X >400 on and v >10^ cm/s. In Is the packet 

® O 

would diffuse to a size of 7x10^ cm, tdiirh would take «^1/4s to convect past 
the antenna system. However, at these dimensions, expanding solitons would 
start interacting with each other. This situation has not been investigated 
and is beyond the scope of the present paper. 

(a) Radiation from Solitons 

A sollton is a spatially confined structure in which intense electrostatic 
fields oscillate at frequencies very close to As pointed out by Papado- 

poulos 8uid Freund (1978), such fields will produce a second order current 
» '*6*’® "•] the density perturbation consistent with Poisson* s 
equation, and X| is the velocity perturbation produced by the sollton electric 
field. If that field has the form (in cylindrical coordinates) 

E.(ilft) * e E(r,z) sin(u t) 

Z 0 

then (Papadopoulos and Freund 1978) 

-e e /E sin(2w t) 

8±rv 

e 

Note that radiation will be produced directly at 2(n^. For the spatial 
structure of a sollton one can take (Degtyarev el al. 1976, Papadopoulos and 
Freund 1978) 
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theory naturally aocounta for the observed predominance of second harmonic 
radiation. 

Perhaps the most convincing evidence of the importance of strong turbu- 
lence effects, and in particular, cf the OTSI, is the ability of the theory to 
account for the observed correlation between the radio and electron fluxes— 
I(2ioe)«Jg*, with aal for weak bursts, and »2.4 for stronger ones. Within the 
context of the theory, bursts which strongly excite the OTSI will amplify 
Langmuir waves well above threshold, into the scaling region where and 

a«2.4; while for weaker bursts and a*1. 

There remains one outstanding problem in the theory of type III bursts 
which we have not addressed; namely, the decay phase. The problem is two- 
fold. In the first place, the decay of the elect ranagnetic radiation is 
observed to be exponential to a high degree and for a long time. The second 
feature is that the product of the decay time and the observed frequency is 
nearly constant over four decades in frequency. Previous attempts to provide 
an explanation for the decay have invoked either collisional damping of the 
plasma waves (Jaeger and Westfold 1949, and Malville 1962), or Landau damping 
(Zaitsev, Mityakov and Rappoport 1972, and Harvey and Aubier 1973). However, 
at decametric and hectometric frequencies collisional damping is slower (by 
nearly an order of magnitude at 100 KHz) than the observed decay rate; and 
were Landau damping the dominant mechanism, radio bursts should cease abruptly 
after the intensity has decayed by only a few e-foldings (Harvey and Aubier 
1973, and Aubier 1974). To develop a complete theory of the decay phase would 
require an understanding of how solitons decay in time, which is an unsolved 
problem in strong plasma turbulence theory, and is beyond the scope of the 
present analysis. 
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FlQilBE CAFTIQMS 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


After Lin et al. (1973)* Electron spectra are shown at various 
times during the type III burst of 16 May 1971. 

Results of a numerical solution of the rate equations that 
describe the OTSI. Parameters were chosen to model the 16 May 
1971 event at 1 AU. For this calculation, L=1.5 AU, nsSxlO”^, 
C=4.6 and f =30 KHz. Panel (a) shows, f„, the distribution 

p€ L 

function of the solar wind and the linearly unstable exciter 
beam, as well as the Langmuir wave energy density, W(+k) — the 
diagonally striped histogram. W(k) has nearly reached W^, the 
OTSI threshold. Panels (b-f) illustrate subsequent stages of 
excitation and stabilization of the OTSI. Ion oscillations are 
depicted by the gray shading. Times are computed from the start 
of the numerical calculations. Calculated values of the 
electron flux are given in panels a, d, and f. 

Similar to Figure 2. Panels (a-c) show results using parameters 

for the 28 February 1972 burst at 1 AU. In this calculation 

C=10.4, L=1.9 AU, n=1.4x10‘‘‘ and f =30 KHz. In panels (d-f) a 

pe 

similar calculation is shown for the 25 May 1972 burst, for 

which c=6.4, L=1.5 AU, n=5.6x10"^ and f =22 KHz. 

pe 

Panels (a-c) show results using parameters for the 16 May 1971 

burst at 0.5 AU. In this calculation L=0.7 AU and f =60 KHz, 

pe 

while n and c were unchanged from the 1 AU calculation. Panels 

(d-f) model the burst at 0.1 AU for which L=0.1 AU and f^„=300 

pe 

KHz. Note that stabilization proceeds as before via the OTSI, 
and reabsorption is unimportant. 
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Figure 5. 


Figure 6. 


Figure 7. 


Figure 8. 


Results at 0.5 AU for two bursts. In panels (a-o) the 25 May 
1972 burst Is modeled using L«0.7 AU and KHs. Panels 

p6 

(d-f) model the 28 February 1972 burst using LsO.9 AU and f^>60 

p6 

KHz. 

Results at 0.1 AU for two bursts. In panels (a-o) the 25 May 

1972 burst Is modeled using LsO.1 AU and f *220 KHz. Panels 

pe 

(d-f) model the 28 February 1972 burst using L*0.2 AU and 
fpgSSOO KHz. Note that stabilization again proceeds via the 
OTSI and reabsorption Is unimportant. 

Yqts equation (III-3.28). Note the 

approximate scaling relations for W-<H(k )<m/M and W(k )>m/M. 
After Fitzenreiter et al. (1976). The electron flux and power 
exponent, a from the relationship I«Jg® are shown for the three 
events for which nvmterical calculations could be performed. 
Observed and computed values of a are plotted. 
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TABLE CAPTIOUS 


TABLE I. Tabulation of tha paranatars uaad In tha nusarical 

Intagration of tha rata aquatlona (X1X-3.4A) - (XXX-3.47). 


T*!a.E I 


Date of 
Burst 

Heliocentric 
Radius of 
Calculation 
(AD) 

Path 

Length 

L 

(AD) 

Plasaa 

Frequency 

f 

I* 

(KHz) 

n- n^/n 

Electron 

^ectral 

Index 

C 

Solar Wnd 
Electron 
Teaperatnre 

(lO^E) 

Figure 

Showing 

Results 

May 16, 1971 

1 

1.5 

30 

5 X lO"® 

4.6 

1.2 

2 (a-f) 


0.5 

0.7 

60 




4 (a-c) 


0.1 

0.1 

300 




4 (d>^) 

February 28, 1972 

1 

1.9 

30 

1.4 X 10"^ 

10.4 

1.2 

3 (a-c) 


0.5 

0.9 

60 




5 (d-f) 


0.1 

0.2 

300 




6 (a-c) 

Hay 25, 1972 

1 

1.5 

22 

5.6 X 10“* 

6.4 

1.2 

3 (d-f) 


0.5 

0.7 

44 




5 (a-c) 


0.1 

0.1 

220 




6 (d-f) 
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